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Abstract 

In this paper, we carry a detailed study of mechanical systems with configuration space 
Q — ► Q/G for which the base Q/G variables are being controlled. The overall system's 
motion is considered to be induced from the base one due to the presence of general non- 
holonomic constraints. It is shown that the solution can be factorizcd into dynamical and 
geometrical parts. Moreover, under favorable kincmatical circumstances, the dynamical part 
admits a further factorization since it can be reconstructed from an intermediate (body) 
momentum solution, yielding a reconstruction phase formula. Finally, we apply this results 
to the study of concrete mechanical systems. 

1 Introduction 

We shall describe a general formalism for studying classical mechanical systems in which some 
of the configuration degrees of freedom are being controlled, meaning that these are known 
functions of time. We will work under the (differential geometric-kinematical) assumption that 
the controlled variables live in the base of a principal fiber bundle Q — ► Q/G = B. The 
remaining variables can be though of as living in a Lie group G and the equations for these 
fiber unknowns are derived by the hypothesis that overall motion respects some (general) non- 
holonomic constraints which are present in the system. 

A special case is that in which the underlying momentum map give conserved quantities, 
even when some of the variables are being acted by control forces. In this case, it is clear that 
motion in the base variables must induce motion in the remaining group variables in order for 
the momentum to be constant during the resultant motion. A concrete example is given by a self 
deforming body for which the shape evolution (base variables) is known and global reorientation 
(group unknown) is induced by total angular momentum conservation (3]. 

In this paper, we consider the more general situation in which fiber motion is induced from 
the base one by the presence of (linear or afHne) non-holonomic constraints. These are repre- 
sented by a distribution D C TQ ([2j Hj) and we shall refer to them as D— constraints. The 
information telling us how the base variables are moving is represented by a base curve c{t) 6 B 
or, equivalently, by a curve do(t) S Q projecting onto c. The desired curve c(t) = g(t) • do(t) £ Q 
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describing the full system 's physical motion is defined by the requirement that it projects onto 
c on the base at each time (i.e. the base variables are the given controlled ones) and that it 
satisfies the corresponding equations of motion plus the D— constraints. The base controlled 
hypothesis can be seen as a set of time dependent constraints and g(t) € G as the do— dependent 
(or gauge dependent, see section 2.4.3| ) fiber unknown. 



The corresponding equations for g(t) are derived by making dynamical assumptions, i.e. as- 
sumptions on the nature of the forces acting on the system. By using variational techniques, 
we give explicitly the equations of motion in section [2] They correspond to the non-holonomic 
momentum equation of [2] with time dependent coefficients evaluated along do(t). Using the 



kinematical structure of the system, in sections 2.4.3|and 4.1 we show how the solution c(t) can 



be factorized by considering specific gauges do(t), yielding that each factor has either a pure 
geometrical (kinematical) definition or it obeys dynamical equations which are simpler than the 
overall fiber ones. 

In section [3j we shall carry out a detailed analysis of systems with a special kinematical 
structure, focusing on the geometric-dynamical factorization of the solution mentioned above. 
Moreover, in section |4j we show that under favorable kinematical circumstances (e.g. in the 
presence of horizontal symmetries [2]), the dynamical factor g{t) of the solution c(t) £ Q admits 
a further factorization. In fact, we can write reconstruction phase formulas |6] for g(t). The 
obtained phase formulas relate the overall system's evolution to the geometry and dynamics 
of simpler intermediate solutions which, in turn, live in smaller spaces (coadjoint orbits). Con- 
sequently, these formulas generalize the ones obtained in [TU] and [3] for rigid bodies and self 
deforming bodies, respectivelly, to the more general setting of D— constrained induced motion. 
Notice that phase formulas become interesting and useful when the dynamical contribution can 
be expressed in terms of the system's dynamical quantities like energy and/or characteristic 
times (see, e.g., [10]) . This is generically accomplished in section [4] and exemplified in section 

El 

The formalism presented in this work, for studying D— constrained, base controlled systems, 
applies to a larger class of mechanical systems than the one encoded in [3J. First, it applies to 
systems with general configuration space Q endowed with a principal bundle structure^ And, 
in the second place, it allows for (linear or affine) D— constraints, and not only momentum 
conservation, to rule the system's dynamics. Indeed, in examples |5.3| and |5.4[ we are able to 
answer two natural questions which arise from [3J: what happens to the corresponding phase 
formulas when magnetic-type forces are acting upon a deforming body, and thus, when the 
(angular) momentum is no longer conserved?; how does a self deforming body move when there 
are additional (internal) non-holonomic constraints between the (no longer controllable) shape 
variables? 

To end this introduction, we would like to comment on the applications of the present work to 
mechanical control theory. First, note that control problems are, in a sense, orthogonal to the 
one we described so far. In this paper, we claim to know the base variables dynamics and we want 
to find the induced fiber motion; while in control theory one starts with a desired fiber dynamics 
and tries to find which base curve induces it (and, after that, how to implement this base motion 
via control forces). Nevertheless, the spirit of this paper is to think of the known base dynamics as 
coming from direct observation or measurement (example |5 . 3| illustrates this point very clearly). 
Indeed, an interesting feature of this kind of systems is the fact that the overall motion c{t) 
can be constructed from that in the base using only the kinematics of c(t) and without actually 



1 Notice that in [3], Q represented specifically the configuration space of a deforming body. 
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knowing the forces which are inducing such base motion. The results on the induced fiber motion, 
obtained by the formalism we describe below, can be thus used for theoretically correcting the 
a priori fiber dynamics prediction when the observed base dynamics deviates from the control- 
theoretical desired one. Also, analytical phase formulas provide interesting tools for directly 
testing different control configurations and theoretical methods. 
Acknowledgements: 
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2 Controlled systems with additional Nonholonomic constraints 
2.1 The Kinematical Setting 

In the remaining, we shall focus on mechanical systems with general non-holonomic linear (or 
affine, see section [33] below) constraints. More precisely, our setting will consists of a mechanical 
system described by the data (Q, L,G, D): 

• Let Q denote the configuration space and G a symmetry Lie group acting on Q by the 
left such that Q — ► Q/G is a principal G— bundle. We shall call, as usual, B := Q/G 
the shape space (see [9]). We denote the action by g ■ q and the induced infinitesimal 
action p gif : TQ — ► TQ. 

• Let k q (-, ■) denote a G— invariant Riemannian metric on Q and k q (-) : T q Q — ► T*Q 
the induced bundle isomorphism. 

• Let L : TQ — ► R denote the G— invariant Lagrangian (with respect to the lifted 
G— action on TQ) given by the (k q —)kinetic energy minus G— invariant potentials (see 
also appendix [6| . 

• Let D C TQ be a constraint distribution. 

We shall assume further: 

(HI) D is G— invariant and D q + V q = T q Q, for all q S Q and Ver q = Ker(ir* q ) denoting the 
vertical subspace of T q Q. This is referred to as the principal case in [2]. 

Now, suppose that, for such a system, the base variables are being controlled in a 
certain known way. This means, that 

(H2) we are given a curve do(t) in Q for t £ I := [ti,i2] or, equivalently, a map c : [ti,t2] — ► 
Q/G s.t. Tr(do(t)) = c(t). The time evolution of the controlled system is then 
described by a curve c(t) G Q such that 

7T( C (t)) = C(t) 

for each t G [ii, t^}. 
The above means that 

c(t) = g(t) ■ d (t) (1) 
for an unknown (do (t)— dependent) curve g(t) in G. 
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Definition: We shall refer to the data (Q, L, G, D, c) as a base-controlled (D— )constrained 
dynamical system. 

The curve g(t) is the unknown for our controlled mechanical problem as stated above. Note 
that, if the controlled problem has a unique solution c(t) for each initial value c(t\) 6 7r _1 (c(ii)), 
then for each curve do(t) in Q lying over c(t) G Q/G, there is a unique satisfying ([I]). In 
this case, the initial condition for the unknown in G reads 

c(t 1 )=g(t 1 )-d (t 1 ). (2) 

The curve do(t) will be called gauge choice or, simply, gauge. This terminology is motivated 
by the analogy between the freedom in choosing among such curves projecting to the same c in 
shape space and gauge freedom in classical gauge field theories (see [9], pj~], references therein 



and also section 2.4.2). 



Remark 2.1 (Restricted configuration space) Note that (H2) implies (but it is not equivalent 
to!) the following holonomic constraint: 

c(t) € n-Hc(I))- 

For a specific problem in which c is fixed, one can restrict the analysis to Q = ir~ l (c(I)). 
Nevertheless, in what follows, we continue with the study of generic c's and thus express the 
results in terms of the kinematical structure of the whole Q. Notice that this is the more 
convenient procedure for studying systems in which c can be (dynamically) perturbed. 

Remark 2.2 (Vertical D— constraints ) Note that the dimension assumption (HI) states that 
the D— constraints are vertical, in the sense that it ensures that the equations of motion (locally) 
drop to the base Q/G with no D— constraints remaining there. In other words, the base curve 
c(t) can be arbitrarily chosen within Q/G. For example, if the sum is direct, i.e., D q ®V q = T q Q 
then D defines a principal connection and we are in the purely kinematic case of [2] . In the case 
Dq Pi V q = but D q © V q ^ T g Q, then constraints are also to be considered in the motion of the 
base variables and, thus, the base dynamics could not be (arbitrarily) controllable. 

2.1.1 Kinematical ingredients 

We now recall some known definitions and properties for mechanical systems that we shall use 
through the paper. 

First recall that for simple mechanical systems with symmetry ([U[Tj) as described above, the 
lifted G action on TQ always has an (equi variant) momentum map J : TQ — ► g* given by 

(J(v q ),X) = (k q (Vq),X Q (q)), 

for X £ q. Let us also recall another ingredients (see ex. |9j): 
• Locked inertia tensor I q : q — ► q* , 

I q = a* o k q o a q 

with a q : g — ► T q Q denoting the infinitesimal generator map, cr q (X) = Xq(c[). Because 
the metric k q is G— invariant, I satisfies the equivariance property: 

I g .q = Ad* oI q o Ad g -i. 
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• The momentum map J is Ad* — equi variant, i.e., 

J(g ■ m) = Ad*J(m) 

with Ad* = (Adg-i) 1 denoting the (left) coadjoint representation of G on q* and * the 
transpose. This follows from the identity 

a g . q (X) = p M (a q (Ad g -iX)). 

Now, from ([TJ we have that 

J f <t) = Pgit^d^d^ig^j^t))) + P g{t ),d {t)j t do{t) (3) 

and thus, 

J(j t c(t)) = Ad^I^g- 1 ^)) + Ad* g(t) J(f t d (t)). (4) 

We can think of Jo(t) := «/(^do(^)) as the apparent or internal momentum along do(t) and 
Io(t) := Id (t) as the locked inertia tensor changing with the gauge motion do(t). 



2.2 Dynamical Hypothesis 

The assumption (H2) above can be interpreted as giving a time dependent type of kinematical 
constraint on the original system, in addition to the one represented by the distribution D C TQ. 
To determine the motion of such a twice kinematically constrained system, i.e. to fmcj^] c(i) in 
Q, we need to add dynamical information. This information consists in assumptions about the 
nature of the forces which are acting upon the system in order to satisfy the imposed kinematical 
constraints. 

For the set of constraints corresponding to the distribution D, we shall assume 

{DH1) D 'alambert 's Principle: The D— constraint forces lie in the annihilator space of 
the kinematical distribution D. 



Denoting F D : TQ — > M. the D— forces (seen as 1— forms on Q) which act on the system 
enforcing the D— constraints, (DH1) means that 

F D (v) = 

for all "virtual displacement" v £ D C TQ. If (DH1) is not satisfied by the system's forces, we 
must then know^jthe D— constraint forces and add them to the equations of motion (see 2.4.1 
below) . 

For the time dependent control like constraints represented by the shape space curve c{t) £ 
Q/G, the assumption takes a less usual form: 

(DH2) The forces which are inducing the motion c{t) to satisfy 7r(c(t)) = c(t) are of a kind 
that we shall denote as good internal ones. Good internal forces seen as 1— forms Ff nt : 
TQ — ► R satisfy 

FLSc) = 

for all vertical variations 5c = Jr| s _ {g{t, s)-do(t)) and any gauge do(t) (see also below). 
2 Equivalently, for a chosen gauge do(t), to find the corresponding g(t) in G. 

3 Or to know some other information about them leading to the corresponding equations of motion. 
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In other words, good internal forces are such that they do not affect dynamically (i.e. by 
adding extra terms) the vertical part of the equations corresponding to (Q, L, G, D). This idea 
is already present in [2] , in terms of the validity of the nonholonomic momentum equations when 
internal (shape space) control forces are present. 

Example 2.3 ( Motion of self deforming bodies) Let Q = M? N ~^ be the configuration space of 
an TV— particles system modeling a deforming body. In this case, usual internal forces between 
the particles of the system satisfying the strong action reaction principle ([5]) are good internal 
forces . For details, see [3J. 

Remark 2.4 (Non good internal forces) If the constraint forces acting on the controlled base 
variables are not of the good internal type, then we must add the extra piece of information 
missing, this is, how the equations have to be modified by adding the non- vanishing terms 



F c {5c) (see 2.4.1 below). In the case of the above example, this means that if there are, say, 
electromagnetic forces acting on the self deforming body which do not satisfy the strong action- 
reaction principle, then one must know the underlying magnetic field data and correct the 



angular momentum conservation equations as usual (see e.g. [5], and also sections 4.5 and 5.4). 



For control purposes, the equations of motion following from (DH1) and {DH2) for the base 
variables r S B = Q/G can be locally written as (for details see |2j) 



M(r)r = -C(r, f) + N(r, f, J(g, g)) + F pot + Fl 



int 



where g denotes the (local) vertical part of variables in Q — Q/G x G, J the (generalized, 
non holonomic) momentum, F pot the potential forces acting on B and F? nt the control forces 
mentioned in {DH2). Also, M denotes the mass matrix of the system, C the Coriolis term 
(quadratic in r) and N a term being quadratic in r and (J, £), where £ is a q dependent element 
in = Lie(G). 

In what follows, we shall assume that the system is being base-controlled, so the control forces 
are inducing via the above equation the prescribed motion c(t) = r(t). The problem is then 
to find the remaining vertical part of the motion, which is induced by the one in the base B 
because of the presence of the D— constraints. 

2.3 The variational principle 

The equations of motion for the above described base controlled D— constrained system, satis- 
fying {HI, 2) and (DH1,2), can be deduced from an adapted variational principle. 
Explicitly, we shall assume that the solution curve c(t) is an extremal of the action functional 

t2 



S Q = J L(c)dt 



ti 

for deformations of the following specific kind: 

c(t,s)=g(s)-c(t). (5) 

These kind of deformations can be called vertical following the ideas of jl]. Also following [3], 
from (DH1) we shall restrict the variations to the ones satisfying the D— constraints, i.e., 

5c e D c( t). 
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Let I = [ii,t2] an d £l(Q;c(t),qi,q2) denote the space of smooth curves I — ► Q with fixed 
end-points qi,q 2 £ Q such that n(c(t)) = c{t). Note that, for a given (any) gauge choice do(t) 
such that n(do(t)) = c(t), any deformation can be written as 

c(t,s)=g(t,s)-d {i). (6) 

and thus 

0(Q; c(t) , q±, q 2 ) « Q do (G; 51 , 52) 

with fid (G; <7i, #2) being the space of smooth curves / — > Q with fixed end points s.t. 
5i • do(U) = qi for i = 1,2. 

So, summing up, our problem is equivalent to the following (gauge invariant) variational 
formulation: 

PI (Gauge invariant formulation) Finding an extremal c(t) of the action Sq, i.e. 5Sq = 
0, among the curves in U(Q; c(i), qi, q 2 ) for vertical deformations Sc(t) = ^| s=0 c(t,s) 

induced by vanishing at the end points, i.e. Sc(ti) = for i = 1,2, and with both c(t) 
and (5c satisfying the D— constraints. 

Once the gauge is fixed, the action S induces an equivalent non autonomous Lagrangian system 
on the G which is, in turn, equivalent to the following: 

P2 (Gauge covariant formulation) Finding extremal curve g(t) in the set f^ (G; g%, 52) for the 
action 

*2 



Sc[do] = J L do (g,g,t)dt 



i.e., for which 55g[^o] = 0, satisfying the gauge induced D— constraints, i.e., 

do(t) + (g- 1 g) Q (d (t)) G D Mt) (7) 
and for variations 5g(t) = ^\ s - g{t} s )) Sg(U) = 0, i = 1,2, satisfying the D— constraints: 



-1 A 
(is 



5) (doC^eD^t). (8) 

s=0 



Note that, although different gauge choices shall lead to different time dependence of the g{t) 
equation 's coefficients, the full solution c(t) is the same for all do 's. In other words, though 
the equations for (and thus g(t) itself) are not gauge invariant, the solution c(t) is. On the 
other hand, g(t) can be seen as being gauge covariant (see remark [2 .5| below). 

In the above formulation (P2), Ld is L(c) with c(t) given by ([!]). It is easy to see that it 
takes the form of the (left) G-invariant non autonomous Lagrangian given by: 

1 ■ ■ 1 / • 

L do (g,g,t) = -k do (d ,d ) + - (l(do(f))£,0 + (J{d )(t),t 

with £ = g -1 <?, I(do(t)) : g — > g* the locked inertia tensor map and J : TQ — ► q* the 



momentum map (recall section 2.1). 
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Finally, note that variations <5£ = ^| s=0 (<7 1 g) induced by variations 5g = ^\ s=Q g{t, s ) 
satisfy the following identity: 

8Z-^{g- 1 5g) = [$,g- 1 6g] (9) 
where [, ] denotes the Lie bracket on q. 

Remark 2.5 (Gauge covariance) Note that (J7| is gauge dependent, meaning that it is different 
for different choices of the gauge curve do(t). Nevertheless, since D is G— invariant, it is gauge 
covariant: if do{t) = g C g{t) ■ do(t) is another gauge, then g(t) in (|TJ satisfies the D— constraint 
equation Q for do(t) iff g(t) := g(t)g~g(t) satisfies the eq. analogous to eq. for the new 
gauge doit). 



2.4 Equations of motion 

Note that, as a consequence of Newton 's second law, the equations for the unknown g(t) shall 
be second order ones. Also, by the time dependent control constraint, they shall also be non 
autonomous and gauge dependent, i.e., its coefficients will depend on time through the chosen 
d (t). 

We shall start with the gauge invariant formulation (PI). Taking into account that the 
variations are of the form (j5j , 

ss Q = o 

straightforwardly implies 

i* c{t) (p(c)) = (10) 

for i c M : C W <—> g denoting the inclusion and 

<*) := {X G g, X Q (c(t)) G D c(t) }. 

The above equation is equivalent to the non-holonomic momentum equation of [2], evaluated on 
the controlled curve c(t) of eq. ([I]) . 



Remark 2.6 (Non necessity of (HI) nor (H2) ) Equation (10) is one of the equations of motion 



of any system whose kinematics is as in section 2.1 without the need of (HI, 2). The only 
dynamical hypothesis needed is (DH1) plus the fact that any other force acting on the system 
(seen as 1— forms on Q) is such that it vanishes when evaluated on vertical variations. What 
these last kinematical hypothesis (HI, 2) add is: that no D— constraints remain on the base 



variables and that these are being controlled, so (10) is the only equation of motion (not of 
constraint) left to solve in the system. 

These are k := dimg c ^ = dimQ c ^ = const, equations coupled to the (dimg—k) number of 
D— constraint equations: 

c(t) G D c{t) . 

Notice that, since the shape space variables are being controlled, that is, since (H2) leave only 



dimQ degrees of freedom, eq. (10) and the above D— constraint equations determine uniquely 
c(t) because of (HI). 
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Below, we shall give more explicit equations for the unknown g(t) by fixing a gauge choice 
do(t) and working in the gauge covariant formulation (-P2)[^J To illustrate on the underlying 
calculation, we shall derive the equations directly from (P2), though they can be also derived 
from (10) using the decomposition ([T|. Let us, thus, evaluate 



= 5S G [d 



o 



I(do{t))S + J(do)(t),6£ 



By eq. ([9| and integration by parts, we have 

*2 



d_ 
dt 



I(do(tm + J(d )(t))+adl(l(do(tm + J(d )(t)),g- 1 5g 



where ad| = — (ac^)* denotes the (left) co-adjoint action. Notice that g 1 5g is arbitrary only 
among variations satisfying the D— constraint Q, i.e. 

g-Hg (t)eg d °^. 



Consequently, for i^ (t) : 



t d (t) 



d_ 
~d~t 



g* denoting the canonical projection, 

I(d (m + J(d )(t)) + ad\ (l(d (m + J(do){t) 



(11) 



must hold. These are k = dimg q (constant Vg G Q) equations of motion for the body velocity 
£(t) = g~ 1 g(t) which are coupled to the (dimg — k) nonholonomic constraint equations eq. 
(|7j) also for 

Before passing to the next section, we give some properties of the subspaces which are involved 
in ( |11[ ) and which follow from the G— invariance of D. Recall that, in general, Q q := {X G 
g, Xq^o) G D q } and i q : Q q ^ g denotes the inclusion. 

Proposition 2.7 The following holds: 



■,9-g 



Ad g g q 



Ad g o i q 



ig-q ° Adg 



Example 2.8 (The purely Kinematical case of [2] ) In this case, D n TOrbc{Q) is trivial and 
thus g q = {0} for all q G Q. Eq. of motion (10) is trivial and the motion of the system is only 



determined by the constraint equation c G D. See also section 4.3 



Example 2.9 (The case of Full Horizontal Symmetries [2] ) In this case, there exists a subgroup 
H C G such that g q is constantly {) = Lie(H) for all q G Q. Then, i^(J(c)) is a conserved quantity 



along the solution c(t). See also section 4.4 



Example 2.10 (The case D = TQ and momentum conservation) When D = TQ and so 
the D— constraints are trivial, equations (10) (equivalently, (11)) imply the conservation of the 
momentum J along the solution c(t). This is the case, for example, of a self deforming body 
which freely rotates around its center of mass with conserved angular momentum (P2 lllj). 



4 We do need (#1,2) for (P2). 



9 



2.4.1 Applied Forces 



In the presence of arbitrary additional external forces F : TQ 
of motion are 

*c(t)(^ J ( c )) = C(t) ° °c(t)( F c(t)) 

where o~ c r t \ '■ — * F c(t)Q denotes the infinitesimal G— action on Q along c(t). 
Now, equation (10) can be rewritten as 

d 



the corresponding equations 



dt 



J{c) = r c (t) 



for a curve T c (t) £ Ker(i*^). This T c (t) is fixed by the D— constraint equations and can 
be interpreted as an external (generalized) torque caused by the forces implementing the 
D— constraints (see example 5.4). 
Within the gauge covariant formulation, the corresponding equations of motion are 

d_ 
dt 



mm + J(do)(t) )+adU I(doV))t + J(do)(t) = ** M t) U*Mt) F <V =: Td o( t )- ( 12 ) 



Assuming that there are no external forces and that (DH1,2) hold, we arrive at eq. (12 ) with 
the forces F = F D representing the D— constraint forces acting on the system. Notice that, 
since (DH1) holds, eq. (11) above follows by projecting via i* do ^ '■ 0* (Q d °^)* ■ If we choose 
a splitting g = g rf o(*) £)d (t) w ^] a pO . g — > gdow ^ e corresponding projector, we get that 
the external torque r^ (t) G Ker i* do ^ present in the r.h.s. of eq. (12) can be also written as" 



P 



do(t) [Pg*d (t) F c(t) 



D 



Expression (12) gives dimg equations coupled to the (diniQ—k) equations of D— constraints. 
Nevertheless, notice that in (12) we have (dimQ—k) new unknowns: the D— constraint forces 
F D . 



2.4.2 Bundle Formulation 

The gauge invariant formulation (PI) and the gauge fixed formulation (P2) of the problem, 
both have as underlying G— bundle Qi — * I which is related to Q — * Q/G by the pull-back 
diagram 

Qi — ► Q 
i i 
I Q/G. 

Moreover, Qj is a trivial G— bundle and the corresponding global sections are the gauge curves 
do(t) projecting to c(t) on shape space. Choosing a section, so Qi I x G, we arrive at the non 
autonomous system on G as described by (P2). 

Remark 2.11 (Relation to 1 — d gauge field theories) The setting above gives a description of 
our time-dependent problem in terms of a 1— dimensional gauge field theory. Here, the fields are 
the sections I — ► Qi I x G and G is the gauge group. See also [9] and references therein. 
Notice that, in this context, the corresponding field theory is not gauge invariant since, actually, 
the problem consists in finding the correct gauge taking do into the desired solution c = g ■ do. 

5 The above expression yields the same Td (t) for any choice of P° since F c r t \ vanishes on D C TQ by (DH1). 
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There is also another set of bundles which are relevant for this problem, specially for the study 
of the equations of motion. These are the vector bundles g D , gf. These are related by the pull 
back diagram 

B? — D 
I I 
I Q 

with g D = U qe QQ q . The vector bundle g D can be also defined as a~ 1 (D), for the vector bundle 
morphism 

a : Qxg — >TQ 
: (q,X)^X Q (q) 

with D C TQ seen as a vector subbundle. Note that bundle gf is also trivial since I is 
contractible. For a given choice of gauge curve do(t), there must exist a smooth curve T(t) £ 
GL(g) such that the set 

{T^X^r/^ (13) 



is a basis of g d °W if {Xj}^™ ( is a basis of the vector space g^ ^ 1 ) C g. 
This is the pull-back (to gf) version of the moving basis formulation of [2j 



Remark 2.12 (Vector Bundle non triviality) From example 2.9 we see that the geometry of 



the bundle g D plays a crucial role in the form of the equations of motion. In other words, the 

I and %) 



geometry of g D enters in the noncommutativity of 4. and i*,.-. in eq. (10). Even though the 



bundle gf is always trivializable, if it is not directly trivial, the need of using time dependent 



sections T(t) enters non-trivially in the equations of motion. See also sections 3.4 and 3.3 where 
this effect is isolated from others. 

2.4.3 Non-Holonomic Gauges 

Recall the constraint equations Q which are coupled to the motion ones (11). Being explicitly 
gauge dependent, a natural question that follows is: Is there a gauge, i.e. a choice of do(t), which 
simplifies these equations? 
For eq. ^ we see that if do(t) satisfies 

^do(f) G D Mt) , Vt G I (14) 
then, (J7| is equivalent to the simpler condition 

m g Q do(t) - (15) 

We shall call a gauge do satisfying ( |14[ ) a non-holonomic gauge and denote it as d^ H . 

Following [2], given the base curve c(t) £ Q/G, a geometrically defined candidate for non- 
holonomic gauge d$ H fulfilling eq. ( |14[ ) is given by the horizontal lift of c(t) with respect to 
the non-holonomic connection. The gauge d$ obtained in this way is defined by 



i* <H I(d$ H (t))i d „ H ) (i* dr J(d» H (t)))=o 
A%f H (d» H (t)) = 
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where i* d ^ H I(d^ H (t))i d N H : Q d °" — ► (fl^")* and A Kin : TQ — ► it denotes a it- valued 1- 



form that projects il q onto itself and has D q as kernel. The subbundle il C TQ can be defined 
to be, at each q £ Q, the (kinetic energy metric) orthogonal complement of (0 9 )q (q) within the 

subspace T q (Orb G (q)): T q (Orb G (q)) = (g«) Q (q)eil q (see [5] for details). 

In this case, the gauge factor d$ (t) of the solution c(t) can be kinematically determined from 
the base-controlled dynamics' c(t). 

In an non-holonomic gauge, we also have 

Proposition 2.13 LetdQ H (t) be a non-holonomic gauge and define the non-holonomic body 
momentum by 

II(t) := /«"(i))(rt(i)) + J(dg B )(JL). (17) 

The following holds: 

• J(c) = v4d*/ t %II(t), so II(t) represents the momentum as seen from the moving reference 
frame defined by g{t) in Q, 

• the constraints read g~ 1 g{t) G rf ° 

• the reconstruction of g(i) from i* NH II(t) is: 

• the equation of motion for II(t) reads 

i* d NB(t) T+ U( V + ad *f \ - 1 ( \ n ^ = °' 

• which is coupled to the constraint equation for II(i): 

i \t)m)-Adn)^ d ° H(t) - 

Remark 2.14 (No constraints and the Mechanical gauge) Note that when D = TQ, i.e., when 
there are no constraints, the non-holonomic connection coincides with the mechanical con- 
nection (see for example |9j and references therein) and, thus, the non-holonomic gauge reduces 
to the mechanical gauge dg /ec/l defined by 

J(dff ech (t)) = 0. (18) 
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3 Special Cases 



3.1 The conserved momentum case 

Here we describe base controlled systems with no additional D— constraints, but whose motion 
is governed by momentum conservation. This case encodes an important class of systems in 
which the fiber motion is induced from the base in order to keep the momentum constant, so 
we shall give a detailed description of the underlying Hamiltonian structure. In section 4.2 we 
shall use this description to the study of reconstruction phases for this systems. 

There are two ways of encoding this conserved momentum case in the general D— constrained 
case described above. One is to think that D = TQ and the momentum J as giving conserved 
quantities due to horizontal symmetries of the whole G (see [2] and sec. |4.4| below). Another, 
is to think 

J(c) = \i = const. 



as an affine constraint on the system (see sec. |3.2| below). Both strategies lead to the same 
results that we shall derive below in a (third possible) direct way, by analyzing the corresponding 
equations of motion. 

Since no D— constraints are present in the system, we only need to assume {H2) and {DH2) 



from sections 2.1 and 2.2 respectively. From these, using the variational techniques of |2.3[ it 



follows that the momentum map J is conserved along the physical motion of the system c(t) £ Q, 
i.e. 

| J(c(t)) = 0, Vt. (19) 

The non autonomous, second order equations of motion for g(t), derived from ( Jl~9| ) 
read 

= adl {t) (I (t)m + Jo(t))+Io(t)j t m) (20) 



where we have denoted 



dV y " v ' dt 



and the initial values (g(ti),g(ti)) must be such that g(t±) -do(ii) = c(ii), ^{g-do)(ti) = c(t\) are 
the initial values of our mechanical problem. We shall now focus on the Hamiltonian structure 
of the equations of motion. 

First, note that I q is a linear isomorphism for each q G Q and defines a symmetric scalar 
product (, ) on g by (X, Y) = (I q X, Y). Let do denote any gauge. If we call 

n = / (f)£ + Jo(t), (21) 

the map sending £ i— ► n, which can be seen as a time dependent Legendre transformation, is 
invertible for all t. In fact, 

C(t) = / - 1 (t)(n- J (t)). (22) 

We also see that 

j{j t c(t)) = Ad; {t) u(t) 
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and, hence, (19) is equivalent to 

d_ 
dt 



n(t) = -ad*i 



n(t). 



(23) 



We will now transform eq. (20) to first order non autonomous equations on T*G making 
use of underlying geometrical structures. Recall that T*G is isomorphic as a vector bundle to 
G x q* via left translations, i.e., by taking body coordinates ([HE]). Also recall the two maps 

L 

Gx g* =tg* given by L(g,U) = Ad*U and ir(g,U) = II. 
We can now state the following 

Proposition 3.1 Let g{t) be a curve in G and II(i) = Io(t)g _1 ^g{t) + Jo{t)- The curve g(t) is 
a solution of (20) iff the curve (g(t),U(t)) is an integral curve of the time dependent vector field 

x( g ,n,t) = (g(i 1 (t)(n- j (t))),-ad*^ {t){u _ Mt)) n) 

onGx q*(~ T*G). In this case, i/L( 5 (ti),n(ti)) = fi, then (g(t),Il(t)) G L -1 (^t) « G /or a// 
t G /. 



Remark 3.2 (Time dependent reduction) Recall that we started with an, a priori, 2 x dimG 
dimensional problem, defined by the non-autonomous second order equation (20) for git). Now, 
due to the conservation of the momentum J, we were also able to reduce the dimensionality to 
dimG = dim(L~ l (//)) because II(t) must be Ad*_ 1 ^fi. See also the next subsection. 



Note that, from the above proposition (equiv. form eq. (23)) we have that 

n(t) g o M c * 

with denoting the G— coadjoint orbit through fi in g*. 
Finally, to solve for g[t) G G, 



1. we have to solve the non autonomous first order differential equation (23) on 
to obtain Il(t) and 

2. then reconstruct g{t) from Il(i) in the G^— bundle L~ l (n) w G — ► O^. 
This last step is studied in section 4.2 below (see also Appendix V?j. 



3.1.1 Hamiltonian structure for the time dependent system 

We shall now add time and energy variables to the above non- autonomous equations on T*G 
in order to get a usual Hamiltonian structure. 

Let us then consider the extended phase space Pe = T*(G x 1) ~ T*G xRxR* with its 
standard symplectic structure S7. By taking body coordinates on T*G, i.e., by trivializing via left 
translations, 

Pg ~ G x g* x I x R*, 
the form £1 becomes, for (g,U,t,E) G G x g* x E x R* , 

n (g,U,t,E) = ~ {d^g^dg) + (II, [g^dg^g^dg]) + (df,dE) . 
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Remark 3.3 (The new variable E) The last term in the above expression, tells us that E is 
the momentum conjugated to time t. Adding this momentum is the usual way of taking into 
Hamiltonian form time dependent systems (see [1]). 

Next, we consider the Hamiltonian function H : P E — ► R, given by 

H E (g,n,t,E) = l - (II, - (nJoHVMt)) + E + T (t) 

with 

T„(t) = K int {t) - \ (l \t)(J (t)), J (t)> 

and where Ki n t denotes the internal kinetic energy defined in Appendix [6] 

The equations of motion corresponding to the Hamiltonian system (Pe,Q l , He) for the 
desired solution curve j(s) = (g,H,t, E)(s) £ Pe are 



1 



g- 1 £g = i Hm-Mt)) 

= ~ a ^oHm-Mt)) U (24) 

ds ' 



f = \ <n, ^(/o-HOJn) - <n, £(i»Ht)Mt))) + £Ut). 



The above third equation, tells us that ^ 



^ and that if we choose as initial value t\(si 



si, 

Equivalently, 



then s = t. Thus, the first two eqs. above become (22) and (23) respectively, 
they say that (g, II)(s = t) is an integral curve of the time dependent vector field on T*G of 
Proposition |3.1| The last equation for E, endows this momentum conjugated to time with a 



physical interpretation in terms of the kinetic energy K(4rc(t)) of the mechanical system on Q 



(see Appendix [6]) : 



E 



-K(-c(t)) + (U(t),I^(t)J (t)). 



Remark 3.4 (Non conservation of Kinetic Energy) Note that though the Hamiltonian H is a 



conserved quantity along the solutions of (24) on P E , it does not represent the kinetic energy of 



the original mechanical system on Q. In general, the time-dependent control forces on the base 
variables do work on the system, implying that the energy is not conserved. Also notice that 
thus, in general, the variable E will not be a conserved quantity, but it will obey non trivial 
dynamics. See also Appendix [6] for more details. 



Remark 3.5 (Mechanical gauge) In the mechanical gauge (18), i.e., when do(t) is horizontal 
with respect to the mechanical connection on Q, we obtain 



E 



Remark 3.6 (Symmetries of (Pe,Q L , He)) The G action on Pe given by 

(g,U,t,E) (hg,U,t,E) 
is Hamiltonian. The corresponding conserved momentum map is 

: (g,U,t,E)^Ad* g U. 
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The corresponding reduced space (|lj) is 

L E 1 (p)/G >1 = 0, 



giving the Hamiltonian configuration for the non autonomous equations (23) on described 
above. 



3.2 Affine D— Constraints 

In this subsection we shall follow [2 a and [3J to show how to handle affine D — constrained con- 
trolled systems. By an affine D— constraint we mean one of the type 

»?(?)= 7(9, *) (25) 

where % D : TQ — ► TQ is a linear fiber projector defining an Eheresmann connection with 
Ker% D = D C TQ. We shall denote, as usual, the vertical subbundle by V = Im% D C TQ. 
The field 7(5, t) is then vertical valued, that is, 7(5, t) E Vg, t. Since our setting involves 
the geometry of the principal G— bundle Q — > Q/G, we assume the following compatibility 
conditions to hold: 

(i) 2l D is G— invariant, that is p gifq o = 21^ o p gJfq , 
(ii) 7 is G— invariant, that is 7(5 • q,t) = p gifq ^(q,t). 

From the G— invariance of % D follows the G— invariance of D. We further assume the dimen- 



sion condition on D, namely, (HI) of section 2.1 Now, we consider the affine version of the 
Lagrange-D 'alambert principle present in [3]: 

PAff The curve q(t) is a solution to the above stated nonholonomic affine constrained system iff 



q(t) satisfies the affine constraints (25 ) and if for any variation q(t, s) with fixed end-points 
such that Sq E D q , then 



5 J L(q,q)dt = 0. 



As in section |2.3[ we adapt this variational formulation to the base controlled case by 
considering only vertical variations c(t, s) = g(t, s) ■ do(t) for some gauge do(t). 

From this, it follows 

Proposition 3.7 The equations for g(t) in order for c(t) = g(t) ■ do(t) to be a solution for the 
affine constrained and controlled system satisfying (i) and (ii) described above are: the equation 
of motion 

ido(t) (| (l(Mm + J(do)(t)) + adl (l(d (tM + J(d )(t)^ = (26) 



with £ = g 1 g, coinciding with eq. (11) for the linear (non affine) constraint case, and the 
constraint equation 



91 D 



do(t) + (g- l g) Q (d (t)) 



l(d (t),t). (27) 
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The fact that the equation of motion for g~^g is the same for the afhne and linear cases is 
already commented, in terms of the nonholonomic momentum equation, in [2] (see page 27). 
As before, we can simplify the constraint equation by choosing suitable gauges oIq. In a 



nonholonomic gauge oIq H , eq. (27) become 



because VI^h {t) (d$ H (t)) = 0. But, if we define an afRne nonholonomic gauge d ^ to be 
one satisfying 

^ {t) {dt ff {t))=l{4 fS {t),t) (28) 



then, eq. (27) reads, 



-1- d Aff (t) 



which is simpler to handle. Notice that eq. (28) plus the requirement ir yd$^ (t)J = c(t) do not 

determine d$ (t) uniquely since dimD can be grater than dimB. On the other hand, when the 
field 7 = 0, a nonholonomic gauge is an affine gauge. 



In section 5.2 we apply this general considerations to study the motion of a controlled ball 
on a rotating turntable. 



3.3 The case G abelian 

We now illustrate on the structure of the equations in the case G is abelian. This allows us 
to isolate the contribution to the motion coming from the non-trivial geometry of the vector 
bundle g D from the Lie algebraic part of the equations of motion (i.e. terms involving ad). 
When G is abelian, Ad g is the identity for all g G G, and thus 

• 5 9 q = 9 £ G, ie, the subspaces fl^are vertically constant in Q, thus Q d °^ = g c W and 

L c(t) - %[ty 

• 1(9 ■ q) = I(q), thus I(c(t)) = I(d (t)), 

• J(c) = I(d (t))g- l g(t) + J(d ) =: II(i), 

• the equation of motion reads 



d 



i Mt)(j t J ( c ))= i 



do(t) 



0. 



(29) 



the constraint equation in a non-holonomic gauge stays as 

g- l 9{t) € = 0^. 



By eq. (16), the constraint equation in terms of J(c) reads 



r\d$ H {t)){J(c)- J(dg u )) = (i* d „ H ol{dS»)oi <B ) (^ r J(c)j e B 



iNH\ 



NH\ 



(30) 
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Remark 3.8 (Base of the g D bundle) Since g 9 ' q = g q for abelian G, the vector bundle g D — > Q 
descends to a vector bundle over the shape space g D — > Q/G. In this context, the objects 
Kit) = K (t) = *£(t) anc ^ -^( c (^)) = I(do(t)) = I(c(t)) really depend on the base curve c(t) £ Q/G. 

Now, we want to re-write the equati on of motion for the momentum J(c) in a usual first order 



differential form. As in section 2.4.2 consider a linear isomorphism Tj : g* — ► g* taking the 
initial fiber ( < H ('O)* to ( «|, 



*< H (t) oT * - Tt ° l *d^ H (t 1 y 



Eq. (29) becomes 

| (»V (t) J(i)) = [ir- 1 ,.^] (j(c)) , (3i) 

which is equivalent to the corresponding expressions in terms of moving basis of |2 . 

The above equation states how the non-triviality of the bundle (g D ) affects the evolution of 

the projected momentum i* dNH ^J{c). Note that even when the bundle is trivializable, but not 
directly trivial, the corresponding equation of motion also contains non-zero TT _1 term. 



Remark 3.9 (Trivial Qj bundle) Recall from section 2.4.2 that the pull back vector bundle 
gf is always trivializable. But, when it is directly trivial, the above equation read 

<l f^H-J(c)) =0 



dt v<%*® 

so it gives a conservation law related to the given base curve c(t). 

1 dtmQ 

gauge curve dQ H (t). Then, constraints (31) for J(c) imply that 



More explicitly, let {e l , NH( .n}*^ 9 ° be a (moving) basis for the fiber g d o (*) along the 



dim S d (t) 



J(c)= mi(d» H (t))e*, +J(d» H ) (32) 



i=l 



for some time dependent coefficients A«(t) G M. to be determined. From (29), we have that the 
Aj(t) 's must satisfy 

A{t) T(t) = -S(i)T(t) - ~c*(i) 
where the time dependent real (dimg d o H ^ x dimg d o H (*)) matrices and -B(i) are defined by 

= ( / ( d o rH ( t )) e dJVff(t)' e ^ff( t )) =: ° 

and the dimg d o H ^ real vector c (i) by 

c i(0 = (^ J ( d v H )> e d* H (t)) ■ 
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Note that A is symmetric and invertible. If we solved these equations for J(c)(t), then the 
reconstruction of g(t) from it is straightforward because, since G is abelian, we can make use of 
the exponential map exp : g — ► G, yielding 

9(t) = exp ^ds lid™)- 1 (j(c)(s) - J(0(*))j (33) 
= exp I J ds (i* d x H o I(d$ H ) o ^v*) (^»( s ) J ( C X S ) 

Vti 

Remark 3.10 (Mechanical connection phase formula) As G is abelian, the above expression 
yields 



c(t) = exp J / ds J(do) (a )J(c)(s) J • gMech{t) ■ d (h) 

with 

9Mech(t) = exp y-J ds Iid^r^id^Xs) 
such that gMech(t) ■ dQ H (t) = HorMech(c)(t) gives the horizontal lift of c(t) £ -B with respect to 



the mechanical connection (18) (see also sec. 4.1). Notice that the equation of motion for J(c) 
(but not the constraint equatiorrl) is the same in any gauge do(t). 



Finally, to better understand how the geometry of the bundle g D enters the equations of 
motion for J(c), we restrict ourselves to the interesting case in which the horizontal space 
with respect to the nonholonomic connection is (kinetic energy metric) orthogonal to the whole 

vertical subspace TOrbc within TQ. In this case, a mechanical gauge do(t), for which J (do) = 0, 



is also a non-holonomic one and eq. (29) yields the parallel transport equation: 



D. ~p=— p i - 7 *V = 0, V l<i<dimg do{t) (34) 

do dt J 

3=1 

for 

ciimg d o(*) 

P*(t):=(j(c),ei, (t) )= £ (WWKw**,® 

3=1 

being the coordinates of J(c) in a basis of g* dual to a basis {e k do ^} dims for which 

(/(d (t))4 (t) ,4 o(t) ) = V 1 < i < dimg d °^, dimg do ^ + !<{< dimg. (35) 



6 For a general non non-holonomic gauge, constraint equation becomes the gauge covariant eq. |7f . 
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Note that, above, for dimQ d °^ + 1 < i< dirriQ then p 1 = by the orthogonality condition (35) 



and because (iff) the constraints (32) are fulfilled. The linear connection coefficients 7I are 
defined by 

dimg 



7fc e d (t)- 



Consequently, for this case, the time evolution of J(c) is geometrically determined, because it 
moves parallel-transported along the base curve c(t) G Q/G in the bundle q d — > Q/G of 



remark 3.8 (see also [2]). 

On the other hand, as noticed in remark |3.9| when the involved geometry is trivial, i.e. 



g D = Q x V with constant V C g, then iyj(c) is a conserved quantity. Indeed, since g is 



abelian, V defines a subalgebra and we are in the case described in section 4.4 



In section 5.1 we apply these general considerations to study the motion of a base controlled 
vertical rotating disk. 



3.4 The trivial bundle case Q = G x B 

To illustrate on how the controlled base variables induce motion on the group variables, we now 
focus on the case in which Q = G x B, i.e., Q — ► Q/G is a trivial principal G— bundle. Recall 
that we are considering the natural left G— action on G x B. In this case, 

TQ = TG © TB 



and thus, by hypothesis (HI) of section 2.1 

D 



'{b,g) 



T h B © 5, 



(M 



with 5, 



T g G n -D(f,, 9 ) as usual. Note that, since D is G— invariant, for each b G B, S(f,,g) 



defines a G— invariant distribution on G which, in turn, is fixed by the subspace S^e) C T e G = g. 
So D is characterized by a smooth map B — > Grdims(&) '■= {Grassmanian of dimS subspaces 
of g} or, equivalently, by a vector bundle 



v = u s, 

beB 



(b,e) 



B. 



(36) 



Conversely, if V — * B is a vector bundle over the base B with fibers 14 C Q, it defines a 
G— invariant distribution D on G x B by setting S^,g) = ^g*eVb- The vector bundle S C D thus 
corresponds to the map 



G x B 

(b,g) 



Gr dimS (Q) 

Lg* e Vb- 



Now, the subspace (recall prop. 2.7) g( b ' 9 ) := {X G g, Xc}(b,g) G is given by 



{XGg, 3Y€S, b , e) ] X = Ad g Y} 



so, 







3 (M 

(6,e) 



g 



(b,e) 



5, 



At this point, we make an assumption on the metric on Q = G x B: 
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(HM) Suppose we have a smooth map 

B — ► {Left invariant metrics on G} ~ {metrics on g} 



O lb ■ 



The metric k®(, ) on Q is assumed to be given by 



for k B (,) being a metric on 5. 

Remark 3.11 (Applicability) This kind of metric on Q = G x B is the one present on typical 
examples (see [2J ) . See also the examples of section [5j 

Assuming (HM), the momentum map J : TQ — ► q* corresponding to the left G— action on 
Q is 

J(b,g) = Ad* g y b (g- 1 g) 

with : q — ► g* denoting the isomorphism defined by the metric (, ) b on q. The inertia tensor 
I(b,g) '■ — > 0* takes the form 



l (b,g) 



Ad* o^ b o Ad g -i. 



Note that we have a natural lift dQ H (t) = (c(i), e) £ B x G for a curve c(t) E £>. This gauge 
dg (t) defines a non-holonomic gauge as defined in section 2.4.3 In fact, this dQ H (t) coincides 
with the horizontal lift of c(t) from (c(ii),e) with respect to the non-holonomic connection of 



[SJ. Moreover, it is also a mechanical gauge (18). 
For this gauge choice, the inclusion 

*< H (t) : d ° W = %(*),e) ^ 
depends only on the base curve c(t) E -B and coincides with the inclusion 

h(t) '■ Vc{t) ^ 



where Vc{t) = S{c{t),e) 1S the fibre of the vector bundle (36). 

The curve c(t) describing the motion on the constrained and controlled system on Q will thus 
be 

c(t) = (c(t),9(t))= 9 (t)-d» H (t) 



and 



J(c) = Ad* g I ( ~ c>e) {g- L g) = Ad* g * m (g- L g). 



In this case, equations of motion (10) read 



d 



J(c) 



or, equivalently, 



d 



*c(t)(0 9)) +ad*_ lg V m (g L g) 



0. 



(37) 
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The constraints for g(t) are 



9 1 9(t) E 



(c(t),e) 



S, 



(g(t),e) 



V; 



£(*)■ 



(38) 



Eqs. (37) can be re- written using a moving basis system on the vector bundle V — * B as 



done in the previous section, yielding the local expression of the non-holonomic momentum eqs. 
of [2] evaluated along c(t). 

Lets simplify the situation a bit more to try to isolate the Lie-algebraic (vertical) contribution 
to the system 's motion from the g D —geometric (horizontal) contribution studied in the previous 
section. 



In case the bundle V 



B is trivial, that is S, 



8(t) 



(b,e) 



l* Vt 



So C 5 for all b £ B, then 



and so equation (37) reads 



d 
dt 



(io*e(t)(5 l 9)) = -i*o (ad* ^ m (g 1 Slj) 



which is an eq. for ^c(t)(9 9)j coupled to the constraint equation (38) for g~ 1 g. Its algebraic 
structure is still hard to handle in general. If we wanted to solve the above (general) equation by 
using usual Lie- algebraic properties of g, then we would need to assume some additional condition 
on how the subspace Sq changes when moving vertically along the fiber (c(t),e) (c(t),g). 

Suppose, then, that So is Ada invariant. It follows that q c ^ = q c( -^ = So and that Sq C q is a 
Lie subalgebra. By the constraints g~ 1 g £ So and eq. of motion (10) becomes the conservation 



law (as in remark 3.9) 



d 



d_ 
dt 



dt 

*o*5(i)(f 



i* J(c)) 



- l g) 



ad* 



_ lg i*o*m(9~ 



l 9) 



Although being integrable, this equation is still hard to solve explicitly in general (see [7] for 
the rigid body q = so(3) case). Nevertheless, in this situation, the dynamical factor g(t) of c(t) 
can be reconstructed from a solution to the above equation in So yielding corresponding phase 
formulas, as described in section [4] and Appendix [7} 

From the analysis of this section, we see that even under very favorable hypothesis on the 
geometry of Q and D, the equations of motion can be very complicated and we cannot continue 
with the general study of c(t). Nevertheless, if we require deeper compatibilities (as above) 
between D and the G— action, e.g. horizontal symmetries, in sees. 4.4 and 4.5 we shall show 



that further phase formulas can be given for characterizing the solution c(t). 



4 Reconstruction and Phases 

In the following, we focus on reconstruction phases (|6]) for both the full solution c(t) and vertical 
(gauge dependent) unknown g(t). The interested reader can find various types of reconstruction 
phases in [8]. 
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4.1 Gauges and phases in Q 



Q/G for D— constrained systems 

c(t 2 ). 



Choice (16) for the non- 



Suppose that the base curve c(t) G Q/G is closed, c{t\) 

holonomic gauge d,Q H {t) provides us with a geometric phase in the motion of the system in 
Q as follows. Being defined as a horizontal lift, dQ H (t 2 ) coincides with the holonomy of the 
associated to the base curve c(t) measured from the initial condition d$ (t\) = c(t\) and with 
respect to the non-holonomic connection. Thus, the corresponding phase formula is 

c(t 2 ) = gDyn(t 2 ) ■ gNH ' dQ H (h) 



with gNH uniquely defined by dQ H (t 2 ) 



gNH 



d,Q (tx) and where gDyn(t) is the solution of 



equations (11) and (15), with = g Dyn gDyn and time dependent coefficients evaluated along 
this gauge d,Q H (t). 

Another geometric phase gup appears when using the mechanical gauge. Let the gauge d,Q H (t) 
be as above and gMech(t) be defined by requiring do{t) := gMech(t) • dQ H (t) to be the horizontal 
lift with respect to the mechnical connection (18) (see [9j) on Q with gMech{ti) = e. If we write 
9Dyn(t) = 9f)(t) • gMech(t), the corresponding equations of motion for the remaining dynamic 
contribution gjj{t) are 



d 

"*(0 \dt 



I(do(t))9l l 9f>) + ad*-! ■ _ I(do(t))g£-g & ) = 



(39) 



the constraint equations ( 15 ) in terms of g^ read 



which are simpler from the original ones (11) because the J(4r(do)) term vanishes by (18). But 



which are more complicated than the original ones for gnyn- 
The relation between the above different gauge phases read 



(40) 



c(t 2 ) = gDynfo) ■ gNH " c{t\) 

= 5d(*2) • 9Mech{h) ■ gNH • c{t\) 
= 9£>(t2) ■ 9MP ■ c(*l). 

with the second geometric phase being gup = gMechfa) ' 9nh- 

Remark 4.1 (Sympifications from different gauges) In the non-holonomic gauge, the constraint 
equations are simpler and, in turn, in the mechanical gauge the equations of motion become 
simpler. One would like to have both simplifications to hold, but this cannot be achieved in 
general since the horizontal lift with respect to the mechanical connection is not horizontal with 
respect to the non-holonomic connection for general D. Finally, we would like to observe that, in 
some situations, we have additional information about the D— constraints and the non-holonomic 
gauge becomes preferable (see, for example, the next sections). 



4.2 Reconstruction Phases for systems with Conserved Momentum 

Now, we shall elaborate on the reconstruction of g(t) for a solution Tl(£) in C g* , as described 
3.1 in case there are no D— constraints. A concrete example of the phase formulas we 



m sec. 



obtain below can be found in [3 for the motion of a self deforming body. 
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Suppose that we have a solution Tl(t) = Ad*_ 1 ^J(c) £ M for eq. (23) with fj, = J(c) = 
const / and that we chose a linear projector P : q -» satisfying 

Ad h oP = Po Ad h . (41) 

From Appendix [7j we know that we can then write 

g(t) = h D (t) ■ g G (t) 

with the geometric phase go being the horizontal lift of n(£) with respect to connection defined 
by P in the G^— bundle G — ► and the dynamic phase hrj 6 G^ defined by 



d 

dt' 



h D h D \t) = p (r^(J(c) - Ad* a j (t)) 



(42) 



with /i_d(£i) = e. The last step follows from eq. (22) for g(t) where I c m denotes the inertia 
tensor evaluated along the physical motion c(i). 

Suppose now that q has an Ad— invariant scalar product (,), as considered in Appendix [7} 



ll*MII 



and {ui}^^ Sfi denote an orthonormal basis with respect to (,) of the vector 



Let u\ 

subspace g M C Q. In this case, equation (42) becomes 



dt 



h D h D \t) 



*(A0ir 

dimgu 



(2K(c(t)) - 2K int {t) + {j {t),I- l {t)J {t)) - (j (t),I-\t)IL(t))) + 



i=2 



(uij^fi) - (uijJ^Adpoit)) 



C (ty 



(43) 



where K represents the kinetic energy of the controlled system in Q (see Appendix [6| . 
When do(t) = d^ ec (t) is the mechanical gauge (18), 



d 
dt 



h D h D \t) 



,(2K(c(t)) - 2K int {t))%n) + Sg*^,/^) 



(44) 



Remark 4.2 (Locked inertia tensor and physical information in ho) The above reconstruction 
phase formula, in the mechanical gauge, relates the dynamical phase ho to the data of the locked 
inertia tensor I c m and the kinetic energy K, both along the physical solution curve c(i) in Q, 
and to the gauge kinetic energy Ki nt {dQ iec ). 

Remark 4.3 (The case J{c) = 0) In this case, the system's motion c{t) coincides with the 
mechanical gauge d^ ech (t) motion because of (18). We thus say that the induced motion c(t) is 
geometrical with respect to the base one c(t) (see also example 4.6 below). 



Remark 4.4 (The case G^ abelian) In this case, /, 



c(t) - *dM™(£ 



and, thus, the only dynamical 



(i.e. non-kinematical) information needed to evaluate formula (44 ) is the system 's kinetic energy 
evolution K{c{t)). In this case, ho(t) can also be easily integrated by means of the corresponding 
exponential map exp : — > G^. 



24 



4.3 Phases for D— constrained, Purely Kinematical systems 

We recall from [2], 

Definition A constrained system (Q, L, G, D) is said to have purely kinematical (PK) con- 
straints if TQ = Ver © D. 

Since D is G— invariant, it defines a principal connection on Q — ► Q/G. Let A D denote the 
corresponding g— valued 1— form on Q. The constraint equation for c(t) then reads 

A D (c) = 



and the vertical eqs. of motion (11 ) are trivial since Q q = for all q. So we have that 



Proposition 4.5 The motion for a base- controlled system (Q,L,G,D,c) for which D defines 
purely kinematical PK constraints, is of geometric nature with respect to c. In other words, the 
solution c{t) is given by the horizontal lift of the base controlled curve c(t) with respect to the 
principal connection on Q — > Q/G defined by the constraint distribution D. 

Corollary If c is closed in [t\, t^, we then have a geometric phase go in the system's dynamics 
associated to the initial value c(t\) and defined by g G = Hol(c): 

c{t 2 ) = go ■ c(ti). 

Example 4.6 (Deforming bodies with zero angular momentum) If we regard J(c) = as a 
D— constraint for the motion of a self deforming body, with J being the angular momentum 
map, then D coincides with the mechanical connection 's horizontal space. From the above 
proposition, we recover the known fact ([II]) that global reorientation g(t) E SO(3) of such a 
body is geometrical with respect to the deformation c{t). 

4.4 Phases for D— constrained systems with Horizontal Symmetries 

We now analyze a geometric-kinematical favorable case leading to phase formulas for the dy- 
namical factor g{t) of c{t). 

Definition ([2]) A constrained system (Q, L,G, D) is said to have (full) horizontal symme- 
tries (HS) if there exists a subgroup H C G such that 

!• £q(<?) G D q yq £ Q when £ G fj := Lie(H) C 0, 

2. (Full condition) S q := D q n T q (Orb G (q)) = T q (Orb H {q)) Vg € Q. 

Condition (2) above states that horizontal symmetries exhaust the whole vertical kinematics. 
The analysis we give below can be extended to the non-full case, i.e. by assuming only (1), but 



we keep hypothesis (2) for simplicity. Example 5.3 below illustrates the non-full case 



For an HS system, the bundle g D is the trivial one Q x h. Since the inclusion map 



f) = Q g becomes independent of the point q, eq. (10) reads 
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Consequently, i^J(c) gives a conserved quantity during the motion of the system as at the end 
This projection itJ(c) can be interpreted as the part of the total momentum map 



3.4 



of section 

which is compatible with the constraints (see also [2]). 

Next, we shall enunciate a few results which follow from the definition of a system with full 
HS. 

Proposition 4.7 The following holds: 

• H C G is a normal subgroup of G, thus f) is G— invariant Ad g \) = f) ; 



if, Ad g = Ad g if,, 



For each q G Q, let Iq 



o I q o iu : f) — > f)* be the restricted inertia tensor, then 
l%. q = Acr g l}Adg-i, VgeH. 



9 1' 



We shall now describe the appearance of phase formulas for the dynamical factor g{t) of 
the motion c(t) of an HS system. First, recall that in a non-holonomic gauge dQ H (t) the 
constraint equation for the body velocity £ = g~ 1 g becomes eq. (15) which, for an HS system, 
reduces to 

g- l g{t) Gf) (46) 

for all t. From the other side, if we consider the non-holonomic body momentum I4(i) £ g* of 

re ha 

rtw = m = (jj)^j (<pi(t) - xv^w)) • (47) 



eq. (17), because of the constraint (46) , we have that 



Thus, eq.(ll), equiv. eq. (45), become 
d 



* W«» = -^-^.(('PW' = - <i b f !(, n (,)-. ! M™(.)))'» n( ' ) ' 

The above expressions are equivalent to 

i;iL(t) = Ad^i^Jic)) 



(48) 



(49) 



because 



i\J{c) = i\Ad* g {U{t)) = Ad* g (ilU(t)), 



by proposition 4.7 The constraint equation (46) can be also put in terms of II(£) as follows 



^(t)(n(t) - JK H (t))) = (iX] - iJJ^W)) e f). 



(50) 



Eqs. (48) and (50), both determine the dynamics of II(t) S g* from the initial value Il(ii) 
J(c) 



/'• 



Now, from (46) and = e it follows that g(t) G if for all t G t 2 ] • Thus, from (49), we 



can deduce that 



ipl(t) G O 



H 

CJ(c) 
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where denotes the H—coadjoint orbit in [j* through the constant element itJ(c). The 

ij* J(c) 

following (commutative) diagram summarizes the relevant geometric situation 
(<MpI(t)) £L-l(0))~fl - 77x(>* 

for the maps L(g,ot) = Ad* a and ir(g,a) = a, (g,a) G H x fj*. Recall that 7 is the mo- 
mentum map corresponding to the left H symplectic action on H x fy* ~ T*H and that 
H ~ L _1 (i? J(c)) — ^ . defines a principal 77„ bundle over the reduced space O u t . , 

as described in Appendix [7J 

Remark 4.8 (Initial conditions) When the initial conditions are g{0) / e in G, so c(0) = 
g(0) ■ dQ H (0), then eq. ( [46| implies that gr(t) = g(0) ■ gii(t) where gij(t) G H is the solution 
corresponding to the initial condition gn(0) = e. Thus, below we shall focus on the g(0) = e 
case. 

We are now in position to apply the usual reconstruction procedure of [6] for the group 
unknown g(t) G H from a solution i*Jl(t) G O h . . Let P : f) — ► f).* ... = Lie I 77„ ... I be a 

i*J(c) \J(c) \ l t ] J\ c ) J 

linear projector s.t. 

P o Adg = Adg o P 

for all g £ H. As described in Appendix |7j P defines a principal connection Ap : T77 — > j( c ) 
and so: 

Proposition 4.9 Keeping the notations introduced above, let n(£) G 0* be a solution of eqs. 



(48), (50) and i^M(t) its projection onto fj* . Then, the corresponding solution g(t) of the recon 



struction eq. (47) which satisfies the constraints (46) luii/i g(0) = e is such that g(t) G 77 Vi G 7 
and 

g(t) = h D (t) g G (t). 

Above, the geometric phase gdt) is the horizontal lift of iZH(t) G H t . from <?g(0) = e with 

ij* J(c) 

respect to the principal connection Ap on the principal H „ —bundle H — ^ H and the 
dynamic phase hp{t) G H^j^ is defined by the equation 

jhDh D \t) = A F (| 5 ) s = p(^)(4)^H^n(t)-^J(^(t)))) (51) 

M°) = e. (53) 

Remark 4.10 (Physical content of hp,) The above dynamical phase hp depends on the (re- 
stricted) inertia tensor 7^ and on the gauge internal momentum Ad*^i^J(dQ H (t)), both as 
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seen from the reference system which is moving along the physical evolution c(i) G Q. Moreover, 
if the non-holonomic gauge choice is the horizontal one (16), then 



dt 



h D h D \t) 



only depends on l\ t y 



Remark 4.11 (The case itJ{c) = 0) In this case, g{t) coincides with the dynamical phase and 
is given by 

rtw = Ai\ H{t) r l i%j{d% H {t)) 



since itll(t) = by (49). Nevertheless, the full motion c(t) is geometric with respect to the 



base one c(t). The reason is that c(t) coincides with the horizontal lift dQ H oi c with respect to 



the non-holonomic connection ([2]) because of equation (16). Notice that this is true for full 



horizontal symmetries, i.e., when the conservation of it J = exhausts the whole vertical eqs. of 
motion (see also [2]). This result generalizes the one of [11] (see ex. 4.6) on the geometric nature 
of base-induced motion for zero momentum systems to the context of D— constrained systems 
with full horizontal symmetries. 

When f) admits an Ad— invariant inner product, the dynamic phase equation can be also 
related to other mechanical magnitudes. 

Proposition 4.12 Keeping the notations introduced above, suppose that is endowed with an 



Ad— invariant inner product (, ) inducing the isomorphism ^ : \f — > f) and let P : t) 

be the orthogonal projector onto fy^j^y Let {m} be an orthonormal basis for fy^j^ 

g(i;j(c)) 
||*(«J(c))|| 



with m = 



. Also, let the non-holonomic gauge d$ H be defined by the horizontal lift (16). Then, 



the corresponding dynamic phase equation becomes 



dt 



h D h- D \t) 



h D (U 



2K(-c(t)) - 2K int (t) 



dira\) 
+ ^ =2 



i;JiC \n i ,(l!> 



c(t)> 



*(*JJ(c)) 
i\J{c))) 



+ 



U; 



In the above expression for the dynamic phase, (Kint) K denotes the (gauge-internal) kinetic 
energy of the controlled system in Q (see Appendix 6 ). As before, -T^ represents the (restricted) 
inertia tensor as seen from the reference system which is moving along the physical evolution 
c{t) E Q . The above formula relates this physical quantities, which are directly involved in the 
system's dynamics, to the phases appearing during the full H— horizontally symmetric motion 
(see Corollary bellow). 

Corollary: Finally, if the solution U(t) G q* is such that ^n(ii) = ^n(i 2 ) then: 

• 5g(^2) is the holonomy of the base path iZH(t) in the H. t bundle H — : > H . with 
respect to the connection defined by P measured from <7g(^i) = e - 
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the solution for the constrained and controlled system c(t) £ Q satisfies the following phase 
relation at time t 2 : 



c(t 2 ) = h D (t 2 )g G (t 2 ) •< H (t 2 ) 



where d.Q H {t 2 ) is the horizontal lift of c(t) with respect to the non-holonomic connection 
([2]), starting from do(t\) = c(ti). 

when, in addition, the base curve c(t) £ Q/G is closed for t £ [t\, t 2 ], so c(ti) = c(t 2 ), then 
d$(t 2 ) = 9g H ' do(t{) where Qq H is the holonomy of the base path c with respect to the 
non-holonomic connection in the bundle Q — > Q/G measured from the initial condition 
do(t\) = c(ti). So, in this case, 



c(t 2 ) = h D (t 2 )g G {t 2 ) ■ g£ H ■ c(h) 



4.5 Phases for systems with dipolar-magnetic-torque type of Affine Con- 
straints 

An interesting special case of affine constrained systems which do not satisfy hypothesis (ii) of 
sec. |3.2| but present reconstruction phase formulas is the following. 



(ii) The affine constraints are of external dipolar-magnetic-torque form, this is, 

s$f «(*)) = i^Acf^u 

for %[ Mech denoting the mechanical connection (see j9] ) . Equivalently, the affine constraint 
can be put in the form 

J(q(t)) = Ad* hM(t) L 

for some given curve hM^t) E G, with /im(£i) = id and the initial momentum value 
L + G 0*. 

The time derivative of the above equation is equivalent to the following non- conservation of 
momentum equation 

±J(q(t)) = adl Mh -iL(q(t)) 

dt n- M h M 

where the right hand side represents a generalized torque of a very special kind. 



In the section 5.4, we shall study the motion of a body with dipolar magnetic moment in an 
external magnetic field which can be described as a system with affine constraints of type (ii) 
above. This justifies our terminology. 

So we now assume (ii) to hold and that we have a base controlled curve c(t). Next, we choose 



the mechanical gauge d$ f ec (i) (18) because D for the above connection form is exactly the 



horizontal space with respect to the mechanical connection. Since constraints represent dimG 
equations, they fully characterizes the dynamics of the group unknown g(t) in c(t) = 5(t)-d^ ec (t). 
Indeed, D defines a principal connection, thus = the zero bundle and so eq. of motion ( |26[ ) 
are trivial, i.e., = 0. These constraint equations in (ii) can be written as 

Ad ^(t)^(t) J <C ec (*) fa -1 ^) = i j ® = const - 

From this, we see that if we call Rhi(t) '■= h~^(t)g(t) G G and TI(£) := I d M ec ^ (g~ l g), then 

Ad R M (t)K(t) = U (54) 
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so n(t) G Ol C 0*, the coadjoint orbit through Lq, for all t. The corresponding equation giving 
the dynamics of II(i) is 

^U(t) = -ad* . n(t) = -ad* / vE(t). 
dt ^ d M e c (i) n(t)-Ad 9 _ 1 / l M/ iM 1 J 

Note that this equation is coupled to the one that defines II(i) from g(t). Nevertheless, recall 
the map 

L : G x g* — * 

L(i? M ,n) = Af^n. 

Equation (54 ) implies that we are in the situation described in AppendixFFland we can thus ap- 
ply the reconstruction procedure ([6]) on the principal G^ — bundle G ~ L (Lq) — > O^ o to ob- 
tain R-M(t) from a solution II(i) G . This yields the phase formula R]\i(t) = R®/ n (t)Rff om (t) 
where the dynamic phase R®f n {t) lies in G^ o and -R^ om (t) is a horizontal lift of n(t) with respect 
to some chosen P— connection Ap in the G^— bundle L~ x (Lq) — G — ► (see Appendix FJ). 
In this case, the dynamic phase equation, when put in terms of the original g(t), reads 

jR D A r i R D h r-\t) = A P (f t R M (t)) g (55) 

= P (RiiWh) - h Mh- h })) (56) 

R D M Vn {h) = ^. (57) 



In section 5.4 we shall work out the details of the above reconstruction formula in the magnetic 
dipole example. 

Finally, if II(ii) = 11^2 ) then we have a phase formula which fully characterizes the motion 
of the system c(t) G Q at time t2- 

c(t 2 ) = h M {t 2 ) • i?f/ n (t 2 ) • Hol^(t l 2 ) ■ gup ■ c(ti) 



where gup is the mechanical-gauge geometric phase (sec. 4.1 ) and Hol^ t ^ ^ is the holonomy of 

the curve II(i) with respect to the P— connection in the bundle L~ 1 {Lq) ~ G — > O^ o measured 
from the initial value e G G. 

5 Examples 

Here we illustrate our general considerations on simple examples of base controlled, D— constrained 
systems. Examples of shape-controlled self deforming bodies with conserved angular momentum 
can be found in [3]. 

5.1 Vertical Rotating disk 

We consider the vertical rotating disk example from [2]. This gives an example of the systems 



considered in section 3.3 In this case, Q = M? x S 1 x S 1 3 q = (x,y,0,ip) and we consider 



G = M 2 x S 1 3 g = (x, y, 9) (left) acting on itself. The Lagrangian reads 

1 .2 .2 1 ■ 2 1.2 

L(x,y,6,<p) = -m(x + y ) + -16 + -Jip 
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and the nonholonomic constraints (non sliding) are given by 

x = RcosipO 
y = RsirupQ 

where R is the radius of the disk. In this case, the base controlled curve is c{t) = ip(t) and 

dg H (t) = (*o, 2/0, 0o, ¥>(*)) 

is a nonholonomic gauge (which, in this example, also coincides with the mechanical gauge). 
Also, 

Q q = span{{Rcos(p, Rsirup, 1) G Lie{G) = M 2 © s(i)} 

and 

J(c) = (mx,my, 10). 
the constraint equation in terms of J(c) for this nonholonomic gauge reads 



3.3 



From section 
I^J(c) G Q", or, 



J(c) = \(t)(mRcosip(t),mRsinp(t), I) 



for some X(t) S 1 to be determined by the corresponding equation of motion (29) for J(c): 



d 



i# do (t) ( fa J ( c ) ) i {Rcos(p, Rsirup, 1 ) 



X(mR 2 + 1) + A 



(mRcostp, mRsinip, I) 



■ (Rcosip, Rsirup, 1) 




0. 



Note that the second term in the last equation is zero because the two vectors are orthogonal. 
Then, since A(i) = by the definition of the momentum J(c), we have 

\(mR 2 + I) = 8(mR 2 + J) = 
which is the vertical equation of motion derived in [2] . The above conservation law can be directly 



computed via equation (34) since 7} = (the underlying linear connection in the 1 —dimensional 
bundle q d — ► Q/G = S 1 is fiat, see sec. 3.3). Consequently, 9 is constant. Finally, since we 



have solved for J(c) using the eq. of motion and of constraints, we can apply formula (33) 
obtaining 



g(t) = I OmR I J ds cosp(s) I ,6mR I J ds cosp(s) ) ,I9(t — t\) 

Note that gMech(t) = (0, 0, 0) in this case. Finally, the full solution c(t) G Q is 
c(i) = g(t)-do(t) 



6mR 



1 1 

J ds cosp>(t 



s) ,BmR\ / ds cosp>(s) ,10{t-t{) ■ (x , y , O , cp(t)) 



t 



OmR I J ds cosp(s) I + xo, OmR I J ds cosp(s) ) + yo, I0(t — t\) + 0$, ip(t) 
.ti I Vti 
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from which we clearly see that motion is induced on the group variables from the base controlled 
curve ip(t) due to the presence of the non-sliding nonholonomic (D— )constraints. 



5.2 Ball on a rotating turntable 

We also recall the setting for describing a ball on a rotating turntable from \2\. This is an 



example of the systems considered in sections 3.2 and 3.4 The corresponding Lagrangian on 
Q = R 2 x SO(3) is 

L = l -m{x 2 + y 2 ) + l -mk 2 {ujl + u? y + w 2 ), 
and the non-sliding affine D— constraints for the ball motion are 

—x + aujy = Sly 
y + acj x = fix 

where (x, y) £ R 2 denote the ball's position and g = u x ^ x (g) +uj y £y (g) +uj z ^f (g) the angular 
velocity of g(t) S SO (3) representing the ball's rotation around its center. Here, ^f(g) denotes 
the right invariant vector in T g SO(3) whose value at e is £j € 50(3), the generator of rotations 
about the i— axis. Also above, a is the ball's radius, mk 2 its (any) principal moment of inertia 
and SI the given angular velocity of the rotating turntable. To take these eqs. to the form of eq. 



(25 ) we define 



4 5 
v„ , c , v„ 



1-5 ^ > „,„,. H , 2 
V \\ \\v \\ 

q\\ 9 



4 5 
v q v q 

l{x,y,g) = Sly - 2 + tlx— — 
\\v„\\ \\v„\\ 



where (,) = (,) R 2 + (>) so (3) denotes the kinetic energy inner product on Q = M? x G with 

d_ 

dy 



G = SO(3) and v\ = -jz + atfig), v 5 = § i +a^{g) in TQ. Note that both D := Ker(A D ^ 



Span{a-^; + ^ (g) ; ~ a '§^ + £x (9) J £f (9)} an d 7 are G— invariant for the natural right action 
of G on Q. Also notice that on the previous sections we considered a left G action on Q, so 
we turn the above natural right action into a left one by defining g ■ (x,y,h) = (x,y, hg^ 1 ) in 
M 2 x G. 

In this case, since shape space B is M 2 , the controlled curve c(t) = (x(t),y(t)) represents the 
position of the contact point between the ball and the table as describing a given trajectory. 
So the problem is to find out how the ball rotates (i.e. to find g(t)) due to the presence of the 
non-sliding affine constraints and to the fact that the contact point is moving in this known 



way (x(t),y(t)). From section 3.2, we know that the corresponding equations for the unknown 



g(t) 6 G are the eqs. of motion (11) and the constraint eqs. (27). Also from that section, 
we know that we can simplify the constraint equation by considering an affine gauge \t) 



satisfying (28). In the present example, gfoJW = Span{Ad q -i^ z } with £ z G 50(3) the generator 



of rotations about the z— axis. Also, the momentum map for the above G— symmetric Lagrangian 
is J(x,y,g) = —mk 2 g~ l g £ 20(3) ~ so*(3). One possible affine gauge choice is 

d$ ff (t) = (x(t),y(t),g Aff (t)) 
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with QAffdAff = \ (~V + + iOk + ^2/)£j/) i- e -> with no z— (spatial) angular velocity 

component. Consequently, the full solution c[t) = ((x(t),y(t), g tot (t)) is written as 

c{t) = g{t) ■ 4 ff {t) = (x(t),y{$, gA f f {t)g-\t)) 

with g{t) satisfying: 

1. (Constraints) g~ l g € g d o / (*) = Span{Ad gAff ^-i^ z } 

2. (Motion) (iJ(x,y,g),Ad {gAff(t)g _ 1[t)r ^ z ) =0 



so(3) 



It is easy to see that, by calling gtot{t) = <7yi//(i)5 _1 (t), eq. (2) above reduces to Jf(c) := 
{dtotdtot^z) so ( 3 ) = const., i.e. the z— component of the (spatial) angular momentum is 
conserved, since 

^j(4,y,d),^e,) = | Ai^^)^ + (j(x,y,5),^- ia d. tot ^e Vso(3) 

= * jf(£) + rofc a (^^,^^-iad. tot9 -,e Vso(3) 

and the second term in the r.h.s. above vanishes. Notice that, although we have a conservation 
law, it is a 1— dimensional one and no non-trivial reconstruction phase formulas for g(t) follow 
from it. 



Remark 5.1 (Conservation due to symmetry) Using remark 2.6, we can easily see, by consid- 
ering G as being only rotations about the z axis and acting by left multiplication on Q, that 
eq. (10) becomes directly the above z— component conservation of the corresponding (spatial) 
angular momentum. Nevertheless, notice that this setting does not give any insight on the 
constraint-base- induced motion g{t). 



Now, from (1) above, we get 
and from (2) that 



9 1 9 = Ad 9- A ) f {tr£* 



uj z = const. 



So, finally, the full base-induced group variable gtot{t) in the full system 's motion c{t) is obtained 
as a product of the two simpler factors gAffifyg -1 ^) described above. 

Note that, in this simple example, the factorization result we obtained following our general 
considerations is the same as what we obtain by proposing the solution gtot(t) = gAf f{t)g~ l {t) 
for the constraints plus conservation eqs. as expressed in ref. [2]: 



gtot = -(-y + Six) £(gtot) + -{x + Sly) £*(gtot) + {const) &(g 



tot) 
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5.3 A non-holonomicaly constrained self-deforming body 

This is an example of a base controlled and I?— constrained system presenting phase formulas 



due to (non full) horizontal symmetries (section 4.4 |2j). The system consists of two rigid spheres 



as in Figure [TJ The small ball is attached to the inside of the big one (holonomic constraint) 
which, in turn, can move freely. The key ingredient is that the first rotates without sliding with 
respect to the second. This last requirement represents a non-holonomic D— constraint on the 
total system and we further assume that no external forces are present. This gives a simplified 
model for a small robot (the small ball) moving inside a space-craft (the big ball). As we shall 
see below, this example generalizes the treatment of [3] by allowing non-holonomic constraints 
to induce total body motion from the arbitrarily controlled (base) variables. 

The configuration space is Q = 50(3) x S 2 x 50(3) 3 (R\,r 2 , R3) defined by requiring 
ri(t) = Ri(t)ri £ IR 3 to be the position of the point i with respect to a reference system 
with axes parallel to those of a chosen inertial one and with origin in the corresponding ball 's 
center (see Figure [l]). We denoted by S 2 the 2— sphere of radius r = ||r2(i)|| = const. In this 
coordinates, the Lagrangian takes the simple kinetic energy form 

L(Ri) = T(Ri) = I (R^R^hRi 1 ^) + \^r 2 (r 2 • r 2 ) + \ (r^ 1 R 3 J 3 R^R 3 ) 

l \ / so(3j A I \ / so(3) 

and the 2 non-sliding non-holonomic D— constraint eqs. (for r 2 o = r z) read 



Adx-i {-RiRi l +R 3 Rz L ),U = Ud 1 [-R 1 R^ l + R 2 R 2 - 1 ),U (58) 

I (Ad .U^ + ^i.d = Ud 1 (- J R lJ Rr 1 + i? 2J R 2 1 ),e 



»( 3 ) v ^2 v - - - - j -yj S0 (3) 

Above, I\ = diag(^m\(r + a) 2 ) , ^3 = diag(^m2a 2 ), with a the small ball's radius, are the inertia 
tensors of the balls with respect to their respective centers in the standard basis of 50(3) 
formed by the generators of i— axis rotations, i = x,y,z, and fj, = r ^ 1 _^ 2 . 

Remark 5.2 (R 2 expressions). It can be easily seen that every expression depending on R 2 
given within this subsection, like the constraints above, is invariant under R 2 {t) R 2 {t)R z {t) 
with R z {t) a rotation about the z— axis. This means that they really depend on r 2 {t) = r R 2 (t)z, 
but we keep the rotational dependence for simplicity. Given r 2 {t) G S 2 , one choice for R 2 (t) is 
given by the horizontal lift in the U z (l) bundle 50(3) — > S 2 (see [3]). 



The distribution D C TQ of tangent vectors satisfying eqs. (58) has dimension dimD = 
dimQ — 2 = 6. 

Now, consider the group G = 50(3) 2 9 (R,gs) (left) acting on Q via 

(R,g 3 ) ■ (Ri,r 2 ,R 3 ) = (RR U Rr 2 , RR 3 g^). 
It is easy to see that both L and D are G— invariant. Shape space B = Q/G can be parameterized 



by elements r 2 ,i 6 S 2 and hypothesis (HI) of sec. 2.1 is satisfied. We also assume (H2) to hold, 
which, in this case, means that the controlled part of the motion is represented by a gauge curve 
do(t) = (e,i"2,i(t),e). If c(t) = (R\(t), R 2 (t), R 3 (t)) £ Q represents the full system's motion, 
then r2i(t) = R^ 1 (t)r 2 (t) represents the position of the CM 2 as seen from a reference system 
with origin at CM\ and with axes rotating with R\, i.e. a system rotating with the big ball. 
Indeed, the full motion can be written as c(t) = (R±(t), R 3 1 (t)Ri(t)) ■ do(t) and note that no 
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z 




Figure 1: The big ball ' s rotation Ri(t) and the position of the center CM2 of the small ball, 
both as seen from refeerence system S, are described by r\(t) = i?x(0 r l0 an d ^(t) = B,2{t)r2o, 
respectivelly. S has its origin at the center CM\ of the big ball and axes parallel to those of an 
inertial frame S. The rotation Rs(t) of the small ball about its center CM2 is described by the 
vector r 3 (t) = i? 3 (*)^30- 
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constraints remain on the controlled variable ri,\ (it can be arbitrarily chosen within B ~ S 2 ). 
Also, notice that from the dimQ = 8 variables, as 2 are being freely controlled, we are left with 
4 equations of motion plus the 2 constraints to solve. 

More physically, the problem is to find the total reorientation of the system R\(t) induced 
by the inside motion. This, in turn, is generated by the inner translational motion do(t) of the 
small ball and followed by its D— induced rotational motion R± (fyR^t), both as seen from a 
system fixed to the big ball, due to the presence of the non-sliding constraints. 

Remark 5.3 (Measurement of t%\) The curve ra i(i) is the one that an astronaut standing in 
the space-craft, modeled by the big ball, would see as the small ball's center moves (see Figure 
[TJ. Consequently, it can also be measured in lab conditions, when the space-craft is attached to 
the floor (and cannot rotate), but when the small ball rehearses the same translational motion 
r 2j i that will occur in space. 

We now turn to the equations of motion. Consider the subgroup H := {(R, e), R £ SO (3)} C 
G. It can be easily checked that F)q = (Lie(H))q C D and that, for q = (R\,r2, R3) G Q, 

Q q = f) Span{Ad R -iAd R2 C 3 z } 

with £2 seen as an element of the second 50(3) copy in Lie(G) = 50(3) ©30(3). The above means 
that we are in the presence of non full [^-horizontal symmetries ([2]). Consequently, 

i\J{c) = hRxR^ 1 + fir 2 (r 2 x r 2 ) x + I3R3R3 1 = hRiR^ 1 + Ad^haR^R? + I3R3R3 1 

in 50(3) ~ s0 *(3) = Lie(H)* is a conserved quantity. Above, A denotes the (Lie algebra) 
isomorphism IR 3 — ► 30(3) and 

h,0 = ^r 2 \ 1 



This horizontal momentum represents the total angular momentum of the system [5]. 

Remark 5.4 (Relevance of the present approach due to constraints) We would like to remark 
that, if we considered only H as symmetry group, as it is done for non-constrained self deforming 



bodies (see [3]), then the D— constraints are no longer vertical (remark 2.2). In other words, the 
corresponding base variables become constrained and it would make no sense to think of them 
as arbitrarily controlled or given. By considering the bigger G instead, we restrict to the smaller 
base variable space which are actually a priori arbitrarily controllable. 



Note that dimg q = 4, so the above conservation law represents only 3 of eqs. of motion (10). 
The remaining equation is 

^ (r^Rs) , Ad iAd R2 £ z ) = 0, (59) 
at v 7 3 7*0(3) 

which tells that there is no angular acceleration of the smaller ball rotation in the CM\ — CM2 
direction. This same effect is observed in the ball on a rotating turn-table example (see [2] and 
the previous section). 



Finally, from section 4.4 we know that we can write (reconstruction) phase formulas for the 
system 's motion due to the horizontal conservations. Below, we summarize the Q— reconstruction 
procedure for obtaining the solution c{t) from the base motion c(i). 
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• We start with dg(t) = (e, r 2 ,i(t), e), and c(t) = (Ri(t),r 2 (t), R 3 (t)) 6 Q representing the 
desired solution. 

• To use the results of the previous sections, we choose a non-holonomic gauge d^ H . We fix 
it by d$ H (t) = (Rl,NH,R^ NH Rl,NH)(t) ■ do(t) with 

[constraints + leg.} - - (-R^nhR^nh + R 3,nhR 3 ^ nh ) = Ad Rl NH R 2 ^R 2 \ 



i^J(dQ H ) = hR3,NHR^NH + ij 1 + ^fii,jvH-R2,i^ 20 ^^flJi J Rl>NHRi i NH+ 

+Ad Rl nh r 2 1 I 2 qR 2 {R 2 ,i = 

with r2,i(i) = rR 2 i(t)z and trivial initial conditions for Ri t NH, i = 1,3. Equivalently, we 
could have chosen the horizontal non-holonomic gauge ( 16 ) leading to the same i^J(dQ H ) = 
equation plus constraint eqs. (|58[) plus one more (involved) equation. 



We now write c(i) = (R,g 3 )(t) • d$ H (t). Notice that, since the horizontal symmetries are 
non-full, eq.(15) for g(t) = (R,g 3 )(t) is non-trivial and yields 

9^93 = A(i) Adjj-^ Ai flliArHjfail ^ 

with A(t) G M to be determined. The corresponding vertical equations of motion for g(t) 
read 



[I)— conservation] itJ(c) 



[Eg. (59)] A 



const = Ad R [I^R-iR - XI 3 Ad RhNH R 2 ^ z ) =: Ad R U\t) (60) 

Ad, 



d_ 
dt 



Ad R -i R 1 R + R 3 NH R,3 NH 



"R^ NH AdR ltNH R, 2 A-, 



S0(3) 



with 



h + Ad Rl NHR21 I 2 flAd 
sponding restricted inertia tensor. 



r-Ra, 



,-i + 1 3 : f) — > f) ~ Pi* ~ so(3) the corre- 



Above, g 3 is s.t. Ri nh R3 ; nh g 3 l = R3,l(t) = Ri 1 {t)R3{t) represents the rotational motion 
of the small ball as seen from a reference system with origin at CM 2 and axes rotating with the 
big ball, i.e., is what an astronaut standing inside the big ball would see (see remark 5.3 and 
l\. Also, A = (g 3 l g 3 ,Ad R -i Ad RlNIlR21 £ z 



Figure 



so(3) 



g 3 g 3 \Ad i R2 ^ 



Mi) 



represents 



a dynamical correction to the (spatial) angular velocity of the small ball in the direccion CM\ — 
CM 2 needed for eq. (59) to be satisfyied from an inertial reference frame. 

Notice that the above equations of motion for R and A are coupled. Nevertheless, in the 
obtained factorization 

c(t) = (R Rl,NH,R R l,NH R2,1,R R3,NH Sj 1 ) 

every element as defined above represents a simpler piece from which the overall motion is 
constructively induced from the known one R 2t i(t) on the base. This shows how we can (geo- 
metrically) take advantage of the kinematical structure of the system for writing the controlled 
solution. Moreover, the global reorientation R can be further factorized by implementing the 
phase formulas corresponding to the f) —conservation reconstruction (sec. 4.4). 



37 



The phase formula for R. From sec. 4.4 we know that R(t) can be reconstructed from the 
body total angular momentum H^(t) solution on Oj*j^ ~ ^radius=i* j(c) ^ Lie(H) = 50(3) ~ M 3 

within the U(l)— bundle 50(3) — > 5 2 _m., m (see details in |3j). In this case, IF(£) was 



defined in eq. (60) and, from mm via 50(3) ~ 



n*(t) = n*(t) x ( (jj 



n*(t) + M3R1 .NH 



r 



with x standing for the usual vector product in M 3 . This equation coincides with the one 
generically presented in [3] but in a very precise non-holonomic gauge d$ H = (e, 33) • rig^, 
which makes the whole procedure compatible with the D— constraints. Also in this case, this 
equation appears coupled another equation, i.e. that of A, since the horizontal symmetries are 
non-full. 

The phase formula corresponding to the reconstruction of sec. 



4.4 



for it J / 0, reads 



i* J 

R(t) = exp I 9 Dyn {t) *> 




Geom 



with the constant i^J G 50(3). The geometric phase R Geom (t) is the horizontal 
total angular momentum curve if^i) in the £7(1)— bundle 50(3) 



5 2 

radius- 



t of the body 
with respect 



to the connection A g (g) 



99 



-1 l l J 



VIA 



(for details, see [3]). 



so(3) 



The dynamical phase 9 Dyn {t) G U(l) = H^*j is defined by (recall sec. 



4.4) 



e Dyn {t) 



i*J 



ds [2K{j t c{s)) - 2K int {s) + \{s) Wrn 2 d> 



Ad 



R2 1 



it J 



so(i) 



Dyn 



|mi(r + a) 2 + ^m 2 a 2 

where if represents the kinetic energy of the whole Q system given in Appendix ^ and l| = 
-^1 + ^2,0 + ^3- Rotation R 2 is defined by r 2 (t) = rR 2 (t)z, giving the physical motion of CM 2 in 
c(t). Notice the unavoidable (dynamical) A dependance of the dynamical phase formula due to 
the fact that the horizontal symmetries are non-full (also compare to the non-D— constrained 
case of |3]). 

Finally, it is worth noting that, when the solution H^(t) is simple and closed for t £ [£1,^2]) 
then 

R(t 2 ) = exp I {e D y n (t 2 ) + e Geom ) s ' ] 



Ri(h) 



with 6 Dyn (t 2 ) as given above and Q Geom given (mod. 2ir) by minus the (signed) solid angle 




enclosed by in the 2— sphere of radius 



within M 3 ~ 50(3). The above is an example 



of a (D-) generalized self deforming body phase formula, not encoded in [3] . 

Remark 5.5 (Control) The above formulas can be useful for control purposes, this is, when 
you want to find the suitable base curve R 2j i(t) inducing a certain desired global reorientation 
R(t 2 ). 
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Remark 5.6 (The case i^J = 0) In this case, the equation for R is geometrical, meaning that 

it is a horizontal lift equation along d^ H = (e, 53) • d$ H with respect to the fj— mechanical 
connection. Nevertheless, this equation is coupled to that of 53 (equiv. A) which is not of 
geometric nature. Consequently, the complete motion induction from the initial controlled base 
variables c(t) = i?2,i(i) £ B is not entirely geometrical. The cause is that horizontal symmetries 
are non-full (compare with remark 4.11 ) and so they do not exhaust the whole vertical dynamics 
(i.e. because of the additional dynamical eq. (59), see also [2] for similar comments). 



5.4 Deforming body with dipolar magnetic moment in an external magnetic 
field 

Here we describe the motion of a (deforming) body with magnetic moment M £ R 3 in the 
presence of an external magnetic field. This system is modeled as an affine D— constrained and 
controlled system for which momentum is not conserved because of the magnetic applied forces 
and which is, thus, not covered by the analysis of [3]. We shall assume the following hypothesis 
about the magnetic nature of the system to hold: 

• the magnetic moment is proportional to the total angular momentum J, i.e. 

M = 7 J 

where 7 is the giromagnetic ratio ([5]). 

• the interaction with an external magnetic field B is of dipolar type ( [5] ) , this is 

— J = M x B 
dt 

where M x B is the external torque acting on the dipole and x denotes the standard vector 
product in IR 3 . 

• the above holds even when the shape c(t) £ Q (see [Hi [H]) of the underlying body and the 
field B{t) are changing with time. 



From the above assumptions, the equation of motion for the angular momentum J(c) of the 
body is 

j t J(c) = 1 J(c)xB(t). 

If we define the corresponding Larmor frecuency vector ([5]) as u>i(t) := — jB(t) 6 M 3 , then the 
above can be re-expressed as 

J(c) = h M (t)L = Ad* hM(t) L 
where hM(t) 6 50(3) is defined by 

hMhjl(t) = Cji(t) 
hM(h) = Id 
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and &i = '§>~ 1 (uJi) G so(3) for the usual Lie algebra isomorphism ^ : (so(3),[, ]) — ► (M 3 , x). 
Also above, Lq denotes the initial value J{c{t\)) seen as an element of so(3)* through the usual 
isomorphisms. 

The equations for the motion of such a system can be derived from the affine-constrained 
Lagragian system {TQ, C, A D , F) where 

• Q — > Q/G is the configuration space of the underlying deforming body ([SJ [3]) with 
symmetry group G = 50(3) 

• the Lagrangian is given by the kinetic energy contribution C{q) = \k q (q, q), where k q a 
G— invariant metric on TQ induced by the standard M 3 — metric ([9]), 

• A D is the mechanical principal connection 1— form on Q — ► Q/G given by 

A D (q) = I~ l J{q) 

where I q denotes the inertia tensor and J : TQ — ► g* the usual angular momentum map, 

• T : Q — > g is the map given by 

F{q)=I q -\Ad* hM{t) L Q ). 
The affine constraints for the physical curve c(t) become 

A D {c{t)) = r(c(t)). (6i) 

We now continue with the analysis in the controlled case, i.e., we add hypothesis (H2) that 
the base curve c(t) G Q/G, representing the changing body's shape, is given. 

Note that the distribution D corresponding to the mechanical connection A D is transversal 
to the group orbit since it is a principal connection (see details in [9J. Then, results from 
section 3.2 in this particular case, say that eqs. of motion (26) for g(t) are trivial (i.e. = 



because g q = OVg). The only remaining equations for g(t) are the constraint ones (27) which, in 
a mechanical gauge d^ /ec (t) (fl~8]) with dQ Iec {t\) = c(t±), read 



Ad * 9 (t) ig d^(t)(9 1 9) = Ad l M (t)Lo 
</(0) = Id. 



Following section 4.5 we call 

R M (t) = h£(t)g(t) e 50(3) 

and note that 

is a conserved quantity. The passage from g to Rm can be understood as passing to describe 
the system from a new reference frame which is rotating via /im(^) with respect to the original 
(inertial) frame, i.e., with spatial angular velocity uji(t) (see [5] pp. 231). 
The above conservation equation can be turned into the form 

L(R M (t),IL(t)) = L e S * 
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with n(i) := I d Me. c ^ {g~ l g) the body angular momentum and L : G x g* — ► g*, L(Rm,^) = 
Ad* RM H. We are in the situation described in Appendix [7j The rotation Ruit) can be thus 
reconstructed from II(t) within the U{1)— bundle L~ 1 (Lq) ~ 50(3) — > ^l " Note that 
0^ o ~ 5 2 for Lq ^ 0. The above defined II(t) must lie in the coadjoint orbit O^ o and satisfies 



n(t) x * ^L w n(t) - Ad p _iw,(t)J . 



-n(*) = n(t) x * I r-^Rm 

The reconstruction procedure follows the lines of [3] and sec. |4.5| Suppose that the solution n(£) 
describes a closed simple curve on the sphere S 2 = , H(t{) = = ^o; reconstruction 

yields 

R M (t 2 ) = ex P ^{e D y n (t 2 ) + e Geom 

where the geometric phase angle Geom can be shown to be (mod 2tt) minus the signed solid 
angle determined by the closed path H(t) on the sphere (see [3]), and the dynamical phase 
6 Dyn (t) is calculated by 




6 Dyn (t) 



U 



ds 



n ( S )' J d Me C(s) n( S ) 



J(c),ui(s) 



In the above expression, the first term gives 2K — 2Ki n t, where K represents the rotational 
kinetic energy (see Appendix [6| and the second term is the magnetic potential energy of the 
system (see [5] pp. 230). Finally, we have a phase formula for the physical curve c(t) £ Q 



c [t 2 ) = h M (t 2 ) ■ exp{B Dyn (t 2 ) 



Lo 



exp(9 Geom -^) -d% Iec (t 2 ) 



L 



which determines exactly the position of the system for the dynamically defined time t 2 in which 
the body angular momentum TI(£) returns to its initial value. This is the affine- constrained 
(magnetic) version of the result obtained in [3]. 



6 Appendix: Kinetic Energy 

Here we derive an expression for the kinetic energy of the mechanical system on Q, in terms of 
the controlled variables curve do(t) and the group unknwon g{t). 

We shall assume that Q is a Riemannian manifold and that the kinetic energy of the underlying 
simple mechanical system (with or without controls) is given by the corresponding metric on Q. 
This means that, if k : TQ ® TQ — ► TQ denotes this metric, then the kinetic energy reads: 

K : TQ — ► R 

Now, on our controlled system, the physical curve c(t) £ Q is of form ([!]) and then, the velocity 
c(t) is given by Thus, the kinetic energy on the controlled curve becomes 

Km) = K int (t) + \ (io(t)m),m) + (Mt),m) 
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where 



Kint(t) = 2 & do(t) (^)(*)>^o(*)) 



shall be called the internal (or gauge) kinetic energy and £(t) = 5 1 ^5'(0 £ 0, io(£) = 
Jo(t) := J(do(t)) as in section 



3.1 



In terms of the body momentum H(t) defined by equation (21 ), the expression takes the form: 



K(j t c(t)) = K mt (t) + \ (l^(t)(n(t)),U(t)) - \ (l^(t)(J (t)),J (t)) 

where the last term can be interpreted as a gauge dependent energy contribution which appears 
because of the use of the "moving reference system" represented by do(t) in Q. 

Remark 6.1 (Mechanical energy) If there are also potential forces present in the mechanical 
system on Q, represented by a potential V : Q — ► M, then the total mechanical energy is 
E = K(f t c(t)) + V{g ■ do(t)). If V is G-invariant, then E = K{f t c{t)) + V{d {t)). Notice that 
as, in general, the control forces are non potential and time-dependent, they do work on the 
system. So the above mechanical energy is not conserved during the controlled motion. 

Remark 6.2 (Gauge potential interaction) In terms of (1-d) gauge field theories, the term 
{Jo(t), £(*)) can be seen as a coupling between the gauge field Jq and the gauge variables £ (see 



also remark 2.11). 



Recall the mechanical connection on Q Q/G (see [9]). The gauge curve do(t) is 
horizontal with respect to this connection iff Jo(t) = for all t. In this mechanical gauge, the 
kinetic energy is given by two uncoupled contributions: 

K{j t c(t)) = K mt {t) + \ <I (t)(«t)U(t)> = K mt (t) + \ (l \t)(n(t)),U(t)) ■ 



7 Appendix: Reconstruction on G — ► M 

Consider the two maps (0) T*G Body ^ oord - Gxq* g* given by L(g, II) = Ad* g Il and n{g, II) = 

II, and suppose that we have a curve (g(t),IL(t)) £ G x q* satisfying L(g(t),H(t)) = fj, = const. 
The idea of this appendix is to reconstruct g{t) from n(£) by means of the fact that n(£) = 
Ad*_!^fi. Note that II(t) lies in the coadjoint orbit C 0* through u. For reconstruction 

([!]), w e need to consider a principal connection on the G^— principal bundle G — > O^, where 
Gfj, := {g £ G; Ad* (J, = fi} denotes the stabilizer subgroup. Recall that this bundle corresponds 
to the reduction L _1 (/i) ~ G — ► O^, where the G^ action on L _1 (/i) C G x 0* is the one 
induced by usual left action (in body coordinates) of G on T*G. Using the principal bundle 
isomorphism 

: L-\u)^UG 
■ {g,Ad* g - 1 u)\ — >g 

we see that a principal connection on G — * can be defined by a choice of a complement 
HOR e C to the isotropy Lie algebra M = Lie(G^), i.e., = HOR e © 0^, and by then right 
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translating this complement to any point g £ G. There is no canonical way of choosing HOR e 
in general. So, let P : g — ► g^ be a linear projector onto g M such that 

Ad h oP = Po Ad h (62) 

for all h £ and define HOR e = Ker(P). The corresponding connection 1-form Ap : TG — ► 
Qfj, induced by P is then given by 

Ap(v g ) g := P{v g g- X ) 

for v g 6 TgG and v g g~~ l denoting the derivative at g of the right translation by g~ l in G. 

Example 7.1 f Ad-invariant metrics) If the Lie algebra q is equipped with an Ad— invariant 
scalar product (, ), then let P to be the orthogonal projector with respect to (, ) onto g M . It can 
easily be seen that this projector P satisfies (62), inducing a principal connection on G — > 0„. 

Now, we shall make use of this connection to reconstruct g(t) from a solution II(t) on the 
coadjoint orbit O^. Following [6]: 

• consider the horizontal lift gc{t) & G from ^g(^i) = g{ti) of the base curve II(t) G 
with respect to the connection Ap, 

• find ho(t) as the curve in G^ fixed by requiring that 

g(t) = h D (t) ■ g G (t) 

be a solution of the reconstruction equation n(£) = Ad*^ 1 ^fi, for the initial value g(t\). 

The group elements in the above decomposition of g(t) at time t, hp>(t) and gc{t)i are usually 
called the dynamic phase and the geometric phase, respectively. The curve hp>(t) must be 
a solution of 

jh D h D \t) = A P ^ t g) g (63) 

with hp>(ti) = e. 

Suppose now that has an Ad— invariant scalar product (,) as in example 7.1 This bilinear 
form induces a vector space isomorphism : g* — > g which transforms the coadjoint action 
into the adjoint action of G. Let u\ = it|^jh an d {ui}f^ Sf * denote an orthonormal basis with 
respect to (, ) of the vector subspace g M C g. Note that this can always be done since ^(/i) 6 g^. 
The orthogonal projector, in this case, can be written as 



and, thus, 



P(v) = Y^ 5ti {ui,v) Ui 



jh D h D \t) =E2*(u i ,^- 1 ) Ui . (64) 
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